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Abstract. COMPTEL imaging analysis revealed a 
patchy, asymmetric distribution of diffuse 1.8 MeV emis- 
sion along the Galactic plane, which is attributed to the 
decay of radioactive 26 A1 in the ISM. If massive stars were 
the major source of Galactic 26 A1, the 1.8 MeV emission 
should be asymmetric and trace the spiral arms of the 
Galaxy, presumed site of massive star formation. Using 
model fits, we indeed find weak evidence in the COMP- 
TEL data that the observed 1.8 MeV emission is at least 
partly confined to spiral arms. We derive a total Galactic 
26 A1 mass of 2.5 M Q from which at least 0.7 M can be 
attributed to massive stars. 
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occurs predominantly inside the spiral arms, especially in 
the case of massive stars. Thus the 1.8 MeV emission pro- 
file should reflect the structure of these arms which is 
thought to be asymmetric with respect to the Galactic 
centre. Previous analysis of COMPTEL data indeed re- 
vealed an asymmetry with more 1.8 MeV emission from 
the southern (Galactic longitude Z=180°-360°) than from 
the northern (Z=0°-180°) Galaxy ( piehl et al. 1995a| ). Ad- 
ditionally, the 1.8 MeV sky map shows lumpy emission 
and 'hot spots'. Prantzos (1993) noted that some emis- 
sion maxima coincide with the assumed tangential direc- 
tions of Galactic spiral arms. Thus a detailed study of the 
spiral arm hypothesis is of interest. In this paper we will 
report on a comparison of COMPTEL phase I+II data 
(May 1991 - August 1993) to models of Galactic 26 A1 dis- 
tribution with special emphasis on spiral structure. 



1. Introduction 

Since the discovery of the 1.8 MeV gamma-ray line emis- 
sion from radioactive 26 A1 by Mahoney et al. (1982), the 
questions of its origin and its distribution along the Galac- 
tic plane stimulated a wave of research (see review of 
Prantzos fc Diehl 1995 ). Core collapse supernovae (SNe), 
Wolf-Rayet (WR) stars, asymptotic giant-branch (AGB) 
stars, and O-Ne-Mg novae were suggested as possible sites 
of significant 26 A1 creation. Early works assumed that the 
large mean lifetime of T2Q ~ 10 6 yr and the low 26 A1 yield 
per source will lead to a smooth and symmetric distribu- 
tion of 1.8 MeV emission. Prantzos (1991) was the first 
who dropped the assumption of an axisymmetric source 
distribution in the Galactic plane if massive stars were the 
dominant 26 A1 producers. He argued that star formation 
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2. Instrument and Data Analysis 

COMPTEL has an energy resolution of - 8% (FWHM) 
at 1.8 MeV and an angular resolution of 3.8° (FWHM) 
within a wide field of view of about 1 steradian. 7-ray 
photons are measured by their consecutive interactions in 
two parallel detector planes where an incident photon is 
first Compton scattered in the upper layer and then ab- 
sorbed (although often not completely) in the lower layer. 
A detailed description of the instrument can be found in 
Schonfelder et al. (1993). We analyzed the COMPTEL 
data in the three-dimensional imaging data-space which 
is spanned by the scatter direction (x,V-') an d the Comp- 
ton scatter angle <p of the incident photons. For the 26 A1 
study, we applied a 200 keV wide energy window, centred 
on 1.8 MeV, to the data. This encloses 48% of all detected 
events from a celestial 1.809 MeV source. 

The Galactic models are represented as 26 A1 source 
density functions. 1.8 MeV model intensity maps are eval- 
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uated by integration of these functions for a grid of Galac- 
tic longitudes I and latitudes b along the corresponding 



lines of sight (Prantzos & Dichl 1995). The convolution of 



these maps with the instrumental point-spread-function 
leads to model distributions of 1.8 MeV source events 
in the data-space. By maximization of the overall data- 
space likelihood, the source models are fitted along with a 
model for the instrumental background to the data. The 
background model was derived using independent mea- 
surements at adjacent energy bands. Data-space analy- 
sis in this approach suppresses continuum emission and 
r eveals only the source s of pure 1.8 MeV line emission 
( Knodlseder et al. 1996 ). To reduce systematic uncertain- 
ties, the unknown (^-profile of the background model was 
adjusted by the fit. 

We tried to eliminate impacts of possible local 1.8 
MeV foreground emission from our Galaxy-wide study: 
for the observed 1.8 MeV emission in Vela, proba- 



bly associated with the Vela SNR (Oberlack et al. 1994, 
Dichl et al. 1995b), and in Cygnus, probably related to 
some nearby star forming regions Qdcl Rio et al. 1994 ), 
two source components with free intensity were included in 
the background model. Both were chosen to have uniform 
intensity within a circle of 10° in radius, centred on (I, 
6) = (259°, 0°) for Vela and on (I, 6)=(83°, 0°) for Cygnus. 
Also, the region of the outer Galaxy between 1—120° and 
/=240° was excluded from the analysis, because the 1.8 
MeV sky-map shows significant emission near the anticen- 
tre which is probably due to nearby 26 A1 sources implied 
by its wide latitude extension. 

The maximum likelihood technique (de Boer et al. 
1992) was applied to determine the parameters and the 
significance of the source models. It uses the parameter 
—2 In A to quantify the model-data agreement, where A is 
the maximum likelihood ratio L (background) / L (source + 
background). Higher —2 In A values signal a better fit of 
the model to the data. Formally, —2 In A obeys a x 2 prob- 
ability distribution, where n is the number of free param- 
eters of the source model. Roughly, the significance of the 
source model over background follows \J — 2 In A, which is 
exact for n = 1. We visualized the model-data agreement 
for each fit by a longitude scan of the data-space for mea- 
surement and fitted model using the software collimation 
technique (Dichl et al. 1993). An acceptance circle of 3° 



was selected which implies an effective angular resolution 
of 10°-12° for the scans. 



3. Axisymmetric models 

3.1. Tracers of Galactic 26 Al 

We start our discussion with axisymmetric models for 
which the source density function only depends on the 
galactocentric radius R and the distance z from the Galac- 
tic plane. For the vertical profile we assumed through- 
out a uniform exponential law with scale height zq. Un- 



fortunately, the Galactic distribution of all 26 A1 candi- 
date sources is poorly known because of the combined 
effects of visual obscuration, uncertain sample complete- 
ness, and small-number statistics. Consequently, we model 
their spatial distribution by tracers observed either in our 
Galaxy or in external spiral galaxies. Many such tracers 
have been proposed up to now, and we discuss only a few 
that seem to be the most appropriate (see also Diehl et 
al., this volume). For this purpose we divide the candidate 
sources in two classes: 

Young population: Stars younger than ~ 10 s yr are 
classified as extreme Population I objects. From the 26 A1 
candidate sources, massive AGB stars, SNe, and WR stars 
fall in this group. The Galactic distribution of these ob- 
jects has a small scale height (zq ~ 90 pc) and can 
be traced either by giant H n regions or giant molecular 
clouds (GMCs). The distribution of GMCs is generally 
inferred from radio observations of the CO J = 1 — > ro- 
tational transition at a wavelength of 2.6 mm. Under the 
simplifying assumption that the 26 A1 emissivity is propor- 
tional to the molecular gas mass surface density we used 
the radial H2 distribution given in Fig. 1 of Dame (1993). 
We added a Galactic centre flat disk component with ra- 
dius of 500 pc to account for the nuclear disk which is 
present in all CO surveys. The total mass of this compo- 
nent was left as free parameter, because the star-formation 
efficiency and therefore the 26 A1 yield in the Galactic cen- 
tre might differ from that in the Galactic disk. 
Intermediate population: Low-mass AGB stars (M< 
2M Q ) and novae have ages greater than a few 10 9 yr. It is 
believed that these objects follow the luminosity profile of 
the Gala xy which is composed of a disk and a bulge com- 
ponent ( Bahcall fc Soneira 1980| ). The disk is assumed 
to be of the exponential form o~(R) cx exp(— R/Rq), 
where o~(R) is the galactocentric surface density and Rq 
is the radial scale length of the disk. Estimates for Rq 
have been derived by numerous workers and span the 
enormous range of 1 to 6 kpc (Wainscoat et al. 1992, 
Kent et al. 1991). Patterson (1984) found a disk scale 



height of ~ 180 pc for the intermediate population. The 
bulge is commonly described as an oblate spheroid with 
an exponential density decrease. Wainscoat et al. (1992) 



fitted IRAS data using p{R, z) cx x 



-1.8, 



where 



\JR? + k 2 z 2 /R c , k c = 1.6 and R c = 2.0 kpc are the axis- 
ratio and effective radius, respectively. However, the rele- 
vance of the bulge component for 26 A1 sources is question- 
able because it possibly contains only objects older than 
10 10 yr, hence no 26 A1 producing AGB stars or O-Ne-Mg 
novae. 



3.2. Results and discussion 

From the fit of an exponential disk we found an opti- 
mum scale length of 5^ kpc and scale height of 180+ 2 gQ 
pc for the Galactic 1.8 MeV emission (all quoted un- 
certainties are statistical 2a errors). The scale length 
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is on the upper side of estimates at other wavelengths. 
The scale height is only weakly constrained because of 
the poor angular resolution of COMPTEL, but it over- 
laps with estimates which reach from 90 pc to 325 pc 
for the youngest and oldest stellar populations, respec- 
tively (Bahcall & Soneira 1980). For comparison of the 



data with the intermediate stellar population, we fixed 
the scale height to 180 pc and added the bulge as addi- 
tional component. The inclusion of the bulge gave only 
marginal, insignificant, fit improvement: there is no evi- 
dence for a bulge on top of an exponential disk in our 1.8 
MeV COMPTEL data (Table |). To test for the young 
stellar population we fitted CO models with and without 
a nuclear disk in the Galactic centre (GC) to the data. 
For both cases the optimum scale height was HO^yQ pc 
which is consistent with the scale height of the molecular 
gas. We find preference for the model with the nuclear disk 
component at the 3<r significance level, but the likelihood 
ratios for both CO models are slightly worse than that of 
the 'best disk'. 

As an illustration for axisymmetric models the longi- 
tude scan of the 'best disk' fit, derived by software colli- 
mation, is shown in Fig. |l| (note, that this scan technique 
has only poor angular resolution (10°-12°), hence sharp 
features in the data are smeared out, large-scale trends 
are emphasized at the cost of not showing small-scale fit 
inadequacies). Obviously, the 'best disk' already gives a 
reasonable first-order description of the 1.8 MeV data. 
However, there are some significant discrepancies between 
the model and the data. While there is a lack of counts in 
the northern Galaxy between / w 50° — 75°, excess counts 
are found around I w —50° and I w —75° in the south- 
ern Galaxy. It is obvious that this north-south asymmetry 
cannot be explained by any axisymmetric model. Actually, 
the 'best disk' is the optimal balance of this asymmetry 
which explains why it yields the best likelihood ratio of 
all axisymmetric models. Therefore, we conclude that ob- 
jects with a smooth, symmetric distribution (as expected 
for low-mass AGB stars or novae) cannot be the solely 
source of 26 A1 in the Galaxy. 

4. Spiral models 

We now drop the assumption of an axisymmetric source 
distribution and investigate the hypothesis that the Galac- 
tic 1.8 MeV emission is correlated with spiral arms. Unfor- 
tunately, the spiral structure of the Galaxy is not well es- 
tablished and even the number of spiral arms is under de- 



bate (e.g. Elmegreen 1985). The most reliable large-scale 
picture of the Galactic spiral structure is probably de- 
duced from the distribution of giant Hn regions. Their 
distances can be estimated spectro-photometrically from 
the distances to the exciting O /B-stars and do not depend 
on models of Galactic rotation (in contrast to spiral struc- 
ture derived from Hi or CO surveys). We adopt the spiral 
model of Taylor & Cordes (1993) who adjusted a four-arm 



Table 1. Fit results for axisymmetric models. The last two 
columns contain the total Galactic 26 A1 mass in the disk and 
in a possible Galactic centre (GC) component. We quote 2a 
statistical errors for the mass which generally increase with 
the number of free model parameters. 





Ro 


Zo 




Disk 
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(kpc) 
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—2 In A 


26 A1 mass (M Q ) 


Best disk 


5* 


180* 


418.0 


3.2 ± 1.1 




-fbulge 


6* 


180 


418.1 


3.3 ±0.5 


0.07 ±0.19 


CO 




110* 


403.8 


2.4 ±0.5 




CO+GC 




110* 


412.3 


2.2 ± 0.5 


0.17 ±0.12 


* parameters optimized by the fit 
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Fig. 1. Longitude scan derived by software collimation for the 
'best disk' fit. The upper panel shows the background sub- 
tracted profile, in the lower panel the residual counts are plot- 
ted (observed-predicted). Note, that we directly compare mea- 
sured to predicted counts, thus an axisymmetric distribution 
must not lead to a symmetric longitude profile because of ex- 
posure variations along the Galactic plane. 



spiral pattern based on giant H n regions and radio-survey 
tangent points to pulsar dispersion measures and inter- 
stellar scattering measurements. Their aim was to obtain 
a quantitative model for the distribution of free electrons 
in the Galaxy to estimate pulsar distances from dispersion 
measures. Recently, Chen et al. (1995) pointed out that 
free electrons could be a valuable tracer of 26 A1 because 
they are mainly produced by the massive star population. 
Therefore we directly compare the Taylor & Cordes (TC) 
free-electron model to the COMPTEL 1.8 MeV data. We 
tentatively added a nuclear disk component (c.f. section 
3.1) of 90 pc scale height in the Galactic centre region 
where the TC model is only weakly constrained by the 
pulsar data. 

The likelihood ratio for both the TC model with and 
without the nuclear disk component is better than that of 
the best exponential disk (see Table |J). From the longitude 
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scan (Fig. ||) we see that the TC model fits the north-south 
asymmetry better than the axisymmetric models. Besides 
the global asymmetry, the TC model also explains the 
two excesses at I w —50° and I w —75° due to the pres- 
ence of spiral arm tangent points in these directions. Ad- 
ditional, but less outstanding tangent points of the model 
at I w —30°, I 30°, and I « 50° are almost invisible in 
the scan because they are smeared out due to the poor 
angular resolution of the software collimation technique 
(see above). However, except the I ~ —30° tangent, all 
spiral arm tangent points of the model coincide well with 
count excesses in the data supporting the hypothesis that 
Galactic 26 A1 is at least partly confined to spiral arms. 

We also studied a more analytical model which con- 
sists of the TC spiral pattern on top of an exponential 
disk. This model allows that some massive stars lie out- 
side the spiral arms, but would also be valid if Galactic 
26 A1 has a composite low- and high-mass star origin. From 
the fit we obtained an optimum disk scale length of 3.5^ 5 
kpc and scale height of 180l~ 2 3Q pc. The likelihood ratio 
is similar to that of the TC+nuclear disk model, thus the 
data cannot tell us which of the two models is a more reli- 
able representation of the Galactic 26 A1 distribution. The 
total 26 A1 mass of 2.7 M is comprised of 2.0 ± 0.7 M 
for the disk and 0.7 ± 0.3 M for the spiral arms. There- 
fore, the total Galactic 26 A1 mass created by massive stars 
is bracketed by ~ 0.7 M Q from the arm component of 
the disk+arm model and by ~ 2.5 M Q from the TC free 
electron model. 



Table 2. Fit results for spiral models. 







Total Galactic 


Nuclear disk 


Model 


—2 In A 


26 A1 mass 


(M ) 


TC 


419.6 


2.6 ±0.3 




TC+nuclear disk 


423.2 


2.5 ±0.4 


0.11 ±0.06 


Best disk+arms 


423.7 


2.7 ±0.8 





5. Conclusions 

We have compared the COMPTEL 1.8 MeV data from 
observation phases I±II to axisymmetric and spiral-arm 
models of Galactic 26 A1 distribution. All models were de- 
tected at a significance level of > 20c above background. 
To first order, the observed 1.8 MeV emission is well rep- 
resented by axisymmetric models. However, details of the 
data like the north-south asymmetry and some regions 
with significant count excesses are better described by 
models which incorporate the Galactic spiral structure. 
Our best fit model holds a total Galactic 26 A1 mass of 
~ 2.5 Mq from which at least 0.7 M Q are produced by 
massive stars. Thus, massive stars clearly contribute to 
the observed 26 A1 in the Galaxy but we can certainly not 
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Fig. 2. Longitude scan derived by software collimation for the 
Taylor & Cordes spiral model including a nuclear disk compo- 
nent. 



exclude from this work that a large fraction of 26 A1 is pro- 
duced by low-mass AGB stars or novae. 
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